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A scalar model of glueball is considered. The model is based on two scalar fields approximation 
for SU(3) non-Abelian Lagrangian. The approach to approximation makes use of the assumption 
that 2 and 4-points Green’s functions are described in terms of some two scalar fields. The model 
is described via non-perturbative method due to value of coupling constant, which does not permit 
us using of Feynman diagrams and therefore of perturbative methods. Asymptotical behaviour of 
the scalar fields are obtained. Profiles of these fileds calculated for a range of values of a parameter 
of the problem is given. Detailed numerical investigation of corresponding equations describing this 
model is performed. The dependence of the glueball mass vs parameters of scalar fields is shown. 
Comparison of characteristics of glueball obtained in our two-scalar model and predictions of other 
models and experimental data for glueball is performed. 
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I. INTRODUCTION 

Nowadays, one of the main problems of quantum chromodynamics (QCD) is the problem of glueball a controversial 
object whose existence is fully admitted, though details of its model are still under question. Glueball represents a 
model for a new hadron which is created by self-interacting gluonic fields. In a sense, a glueball can be imagined as 
a blob of gluonic fields, and this causes difficulties in working with the model. Particularly, the model of glueball 
can be described only in the realm of QCD, since in this case we deal with gluonic fields. The non-dimensional 
coupling constant of SU(3) nonabelian gauge theory is the base for QCD, and the latter being strongly nonlinear 
theory requires the constant to be greater than 1. Thus, mathematically, this condition stipulates that we can not 
use Feynman diagram method, namely, the perturbative quantum field theory. The matter is that the method can 
be used only in case of weak interactions, in which dimensionless coupling constant < 1. Therefore, the problem lies 
in that for description of the model we must use non-perturbative methods of quantum field theory. For example, in 
Ref’s lILdl the authors consider the glueball decay rates in the Witten-Sakai-Sugimoto model, a holographic top-down 
approach for QCD with chiral quarks in Witten’s holographic model of nonsupersymmetric Yang-Mills theory. 

At the present time, experimental search of lightest glueball is conducted based on assessments of its mass, which 
lies in the range of 1000-1700 MeV (see Ref. (3) for the experimental and theoretical status of glueball). Although, it 
is necessary to mention that experimental search did not reveal any strong evidence of existence of glueball. 

In Ref’s Si an approximated approach was proposed for the non-perturbative quantisation in QCD. In this 
approach given SU(3) Lagrangian is turned into Lagrangian of two scalar fields using some assumptions. These fields 
describe two and four point Green functions entering the initial Lagrangian. Hence, we obtain simplified description 
of gluonic fields in the form of two scalar fields. 


II. SCALAR TOY MODEL OF GLUEBALL 

A glueball is a hypothetical composite particle that consists solely of non-Abelian SU(3) gauge field, without valence 
quarks. The existence of a glueball is consequence of the self-interaction of gluons within quantum chromodynamics. 
Nonlinear self-interaction of gluons in QCD leads to possibility of the existence of a color-neutral state made of gluons 
only, which was called glueball. Glueball is also thought as a bound state of gluons, and it’s properties cannot be 
described within a perturbative approach to QCD. Glueball remains an obscure object over thirty years after QCD 
was used to predict such a state. It is well known that the gluon condensate, from which glueball is thought to be 
made of, can only be determined in a nonperturbative formulation of QCD. So far various attempts have been made 
to determine gluon condensate from first principles 0, Q. We refer the reader to Ref. Q for more details. 
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Our main idea is to write an effective Lagrangian. This Lagrangian is obtained from the Lagrangian of SU(3) 
non-Abelian gauge theory. In order to do this, we first separate SU(3) color degrees of freedom into two parts: 
subgroup SU(2) C SU(3) and coset SU(3)/SU(2). Then we average the SU(3) Lagrangian using some assumptions 
and approximation. Our approximation is based on the main assumption that the 2 and 4-points Green’s functions 
are described in terms of some scalar fields </> and \ due to the following relations: 


(G 2 )^(x,x) = (A-(x)At(x)) 

~rnf ~ C“lcj) 2 (x), 

(1) 

(■ F 2 ); b „(x,x) = (dMx)d„A b 0 (x)) 

~ D a b p d n<t>( x )dv<P(x), 

(2) 

( G 4)“tpt (x,x,x,x) = (A“(x)A b v (x)A c p (x)Ai(x)) 

n(A“(x)A b v (x))(A c p (x)Aj(x)) 

(3) 

(G 2 )™(x,x) = (A™(x)A" v (x)) 

~C™ X 2 (x), 

(4) 

( f 2 )™ (x,x) = (d li A^(x)d v A^(x)) 

nD^d^d^x), 

(5) 

(Gt)™£(x,x,x,x) = (AZ(x)AXx)AV(x)A*(x)) 

~ (A™(x)A™(x) - ml) (A p p (x)A q a (x) - mf) - mf, 

(6) 

(G 4 ); b ™(x,x,x,x) = {Al{x)A b v {x)A™{x)Al{x)) 

« CZ b ™ct> 2 (x)x 2 (x) 

(7) 

where a, b,c,d = 1,2, 3 are SU(2) indices, to, n,p,q = 4, 5,6, 

7,8 are coset indices, and Cfj mn , D ab j mn , and 

mi, 2 


are closure constants. Similarly, we see that for turbulence modeling we have to introduce some closure constants. 
An effective Lagrangian then becomes 


g 2 c _g 2 (c \ 

-C cS --{C su{3) ) 


- Y (</> 2 - mj ) 2 + 


A2 

4 




( 8 ) 


where g is the dimesionless coupling constant; \i, 2 , mi ,2 are some parameters; the signature of the spacetime metrics 
is The effective Lagrangian © is an approximation to the nonperturbatively quantized SU(3) gauge 

theory. 

Quantities entering the Lagrangian ([5]) have the following meanings and origins: 

• the scalar fields (j) and x describe nonperturbatively quantized SU(2) and coset SU(3)/SU(2) degrees of freedom, 
correspondingly; 

• the terms V M 0V M </> and V^x^X are the result of the nonperturbative quantum averaging of {S7 p A B ) 2 in the 
initial SU(3) Lagrangian; 

• the terms </> 4 and \ 4 are the result of the nonperturbative quantum averaging of f ABC f AMN A B A < ^A Mll A Nl '\ 

• the term cj) 2 x 2 is the result of the nonperturbative quantum averaging of f Aab f Amn A p A b A rntl A nl '; 

• mi ,2 appear to be closure coefficients; 

• the terms <jrm\ , x 2 m§ arise due to the closure coefficients mis¬ 
using the Lagrangian |8|), we derive the associated field equations describing glueball in the following form: 

d^cj) = -</> [x 2 + Ai (c/) 2 - to 2 )] , (9) 

d^ X = -x[0 2 + A 2 (x 2 -m 2 )]. (10) 

In Ref. Q the spherically symmetric solution to equations (HI) and m is considered. This ball serves as a nonper¬ 
turbative scalar model of a glueball. It is shown in Ref. [§] that such solutions with Unite energy do exist. 

Here we want to investigate these solutions in more details. 


III. TWO SCALAR MODEL OF GLUEBALL 


Our model is described by a system of equations for coupled scalar fields (l9l)- (fT0l) . In order to describe glueball we 
consider the spherical symmetric case. Let us rewrite the equations in the dimensionless form 


= ^[X 2 + Ai (^ 2 -m 2 )], 

x" + ~x' = X[4> 2 +A 2 (x 2 -ml)], 


( 11 ) 

( 12 ) 
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here redefinitions were made 4> —t 4>/fa 0), x —■► x/0(O)> m i ,2 —> /0(0), and the dimensionless coordinate x = rfa 0) 

is introduced. The field equations (EMUD can not be solved analytically and therefore we will solve them using 
numerical methods. 


A. Behaviour of the scalar fields at the origin and at infinity 

Numerical analysis shows that the asymptotical behaviour of the scalar fields is as follows 

fax) —)■ mi,x(x) -A 0. (13) 

Taking into account the asymptotical behaviour m and equations m and m we can find the asymptotical 
behaviour of the scalar fields 


where 


where 


e ~Xy/2\i m\ 

fax) ft 17li (poo -, (14) 

X 

g-rc^/mf —A 2 m| 

Xfa) ~ Xoo -, (15) 

X 

<fioo, Xoo are some constants. Near the origin, the scalar fields change according to the following laws: 

fax) = fa + ^y- + . . • , (16) 

X(x) = Xo + —h • • • (17) 

fa = fa [Xo + A i {<t>l ~ mf)\ , (18) 

X 2 = Xo [Xo + a 2 (xo - m l)] • (19) 


B. Numerical solution 


The numerical analysis shows us that regular solutions of (ED-d exist for some special case of parameters mi ,2 
only. It means that we have to consider these equations’ set as a non-linear eigenvalue problem for eigenvalues TOi ,2 
and eigenfunctions fax),\{ x )- The solution to the system of equations (fTT1) - ([T2l) was obtained via ’’step by step” 
method, when approximate solution is improved in every following step. 

Firstly, we expanded the system of equations near the origin using Taylor series, and got the boundary conditions 

m = 1, ^'(0) = 0, x(0) = xo, x'(0) = 0. (20) 

The system of equations (fTT1) - (fT2l) was solved as a non-linear eigenvalue problem in which values of the following 
parameters were sought: 7711 , 2 - Other parameters except xo i n the system of equations were set to appropriate values: 

Ai = 0.1, A 2 = 1.0, (21) 


while xo was varied in the range of values from 0.05 to 50. For every xo we got appropriate eigenvalues of mi j2 which 
are given in Table Q] Profiles of solutions for a range of xo values are given in Fig. [l][2] 

The dimensionless energy density of the presented solution is 


9 2 1 


1 


1 


ifafa = e{x) = -ft (x) + -/' (x) + y (<j?{x) - m\) + y/ 2 (x) (/ 2 (x) - 2mf) + -/ 2 (x)^ 2 (x) 


( 22 ) 


Thus the dimensionless glueball energy is 

OO 

9 2 1 


Tic fa 


W = W = ■ 


to j x 2 j/' 2 + + £(*’- ml) 2 + (/ 2 - 2m 2 ) + 


dx 


(23) 
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FIG. 1: Bunch of rf)( x ) curves calculated for a range of 
values of xo from 0.05 to 1. The bottommost curve corre¬ 
sponds to the value of xo = 0.05, every curve higher than 
this corresponds to greater values of xo ending up with 
the topmost curve for xo = 1- 


Xo 

mi 

7712 

0.05 

1.01077 

1.01019 

0.1 

1.03487 

1.03189 

0.2 

1.1 

1.08769 

0.3 

1.17779 

1.1501 

0.4 

1.2612 

1.21803 

0.5 

1.3509 

1.28827 

0.6 

1.44491 

1.36096 

0.7 

1.54239 

1.43576 

0.8 

1.643 

1.51235 

0.9 

1.74629 

1.59053 

1 

1.85201 

1.67014 

2 

3.00371 

2.52084 

5 

6.9365 

5.32339 

10 

14.0095 

10.2078 

15 

21.307 

15.1572 

25 

36.1893 

25.1091 

30 

43.7127 

30.0954 

50 

74.1089 

50.0653 


TABLE I: The eigenvalues mi ,2 calculated from 
equations GD-aa for corresponding values of xo- 



FIG. 2: Bunch of x( x ) curves calculated for the range of 
values of xo from 0.05 to 1. It can be seen from initial 
value of a curve that the topmost curve has been cal¬ 
culated for Xo=l) the bottommost curve corresponds to 
Xo=0.05. 


Xo 

X0,o 

x x,o 

0.05 

12.1963 

17.1703 

0.1 

7.39187 

9.61754 

0.2 

4.85681 

5.78707 

0.3 

3.92591 

4.46546 

0.4 

3.3682 

3.7237 

0.5 

3.01152 

3.2655 

0.6 

2.75011 

2.94117 

0.7 

2.54635 

2.69574 

0.8 

2.38164 

2.50203 

0.9 

2.24421 

2.34373 

1 

2.1271 

2.21116 

2 

1.4691 

1.50183 

3 

1.15977 

1.18304 

5 

0.839859 

0.856706 

10 

0.520063 

0.530226 

15 

0.385755 

0.392831 

25 

0.260804 

0.265031 

30 

0.22602 

0.229494 

50 

0.150125 

0.152079 


TABLE II: Characterstic sizes , £ x ,o of regions 
of 4>(x) and x( x ) fields correspondingly calculated 
for different values of xo- 


The glueball energies were calculated for each set of xo together with eigenvalues given in Table [H and plotted to 

Fig. El 

Also, we have calculated the characteristic sizes £</>, o,x Xt o (they characterize the glueball size) for every value xo in 
the following way 


x( x x,o) 


m + 


mi - 0(0) 

2 


x(Q) 

2 


These values are presented in Tabic IXTl and plotted to Fig. 01 


(24) 

(25) 
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FIG. 3: The dimensional glueball energy calculated for 
values XOi mi ,2 given in Table [TJ 



FIG. 4: Characteristic sizes of regions of (j>(x ), x( x ) fields 
calculated for different values of xo given in Table HH Dot¬ 
ted curve corresponds to the characteristic sizes of (j>(x ), 
while crossed one corresponds to the sizes of x( x ) field. 


IV. COMPARISON OF NUMERICAL VALUES OF OUR MODEL WITH THAT OF OTHER MODELS 

Now we want to compare the energy of glueball of our model with that of obtained in other models. We considered 
the models reviewed in Ref. 0 which gives the mass for the lightest glueball within the range of 1000 - 1700 MeV. 
For comparison, we considered one of them with the mass of 1500 MeV. First, we calculated (j >o parameter from the 
energy taken from Ref 0. After that, we calculated the same constant from characteristic size of a glueball which 
supposedly is the radius of a proton, i.e. 1 fermi, Eq. m- 


m = 

g rn g 

h c e’ 

(26) 

m = 

Xq 

ro ’ 

(27) 


where m g =1500 MeV is the glueball mass; g ~ 1; Xq is dimensionless characteristic radius of a glueball taken from our 
model, which is roughly taken equal to 10, we have taken this value because this number represents the characteristic 
size of the region where scalar fields c/>(x), x( x ) are concentrated (i.e. radius of glueball), see Table ITT1 Fig. S] for 
details; E is dimensionless energy calculated for given eigenvalues m 1,2 and parameter xo; r o is the numerical value 
of radius of proton, namely, 1 fermi. 

The value of 4>o we obtained for the mass of 1500 MeV and for x'o=0-05, <f> 0 ~ 10 15 cm -1 and (f>o calculated from 
the data of our model is 4>o ~ 10 16 cm -1 . We can see that these values are in good agreement, taking into account 
that we compared the values qualitatively. 


V. CONCLUSIONS AND DISCUSSION 

Thus, we have investigated the solutions describing glueball approximately in the approach of two scalar fields. 
We have obtained the set of regular solutions having finite energy. Every solution has a good asymptotic behaviour, 
bringing us to the finite energy of this configuration of </>( x),x{x) fields. We have shown that when xo parameter is 
decreased, the region in which (f>(x),xi x ) fields are concentrated grows and correspondingly when the parameter is 
increased the sizes of the region decrease. Physically, this means that in our approximate model the size of glueball 
depends on Xo parameter in the way we described. We have also investigated the dependence of glueball energy on 
Xo parameter. We have shown that when xo — I 0 the energy decreases and correspondingly when xo -A 00 the energy 
blows up. 

It is significant to note that having dependence on xo parameter of the energy of glueball, one can investigate 
statistical properties of glueball which has thermal contact with thermostat. In such case, glueball will play a role 
of statistical object , in which fluctuations of energy occur due to thermal contact with thermostat. Specificity of 
such investigation is in that we have statistical quantum object consisting not of quantum particles but of fluctuating 
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quantum fields. From mathematical point of view this means that we have to calculate statistical sum for non- 
perturbative quantised object. Such problem represents a complicated issue as well as any other problem in the area 
of non-perturbative quantisation. 
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